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TECHNICAL MEMORANDUM 


ANALYSIS TECHNIQUES FOR RESIDUAL ACCELERATION DATA 

L INTRODUCTION 

There are various aspects of observational data that may be of interest to an 
investigator, e.g., mean, variance, and minimum and maximum values. Observational data 
such as biomedical data, economic data, seismic data, and accelerometer data are recorded 
as either con tin uous time functions or discrete time series. While statistics such as those 
mentioned above can be obtained from data in this form, additional information can often 
be obtained by looking at the data from a different perspective, such as can be obtained by 
transformation of data into a different domain or into different coordinate axes. Of 
particular interest to us is the analysis of residual acceleration data collected in orbiting 
space laboratories. A thorough understanding of such data and the ability to manipulate the 
da ta will allow the characterization of orbiters so that investigators can better understand the 
results of low-gravity experiments. 

Most time functions can be considered as the sum of sinusoidal and co-sinusoidal 
terms of various frequencies (all harmonics of some fundamental frequency) and a steady 
state term. Spectral or Fourier analysis consists of a transformation of data in the time 
domain into the frequency domain in which the relative strengths of the frequency 
components of the data can be studied. Although the equations used to perform such 
transformations were originally applied to continuous, periodic functions, the theory can be 
adapted to the analysis of discrete, aperiodic time series. 

The orientation of recording axes can be an important factor in the collection of 
observational data. It may be of interest to an investigator to look at data in terms of a set 
of axes other than that in which the data were collected. Such analysis is performed by 
means of a transformation of coordinate axes in which data in a given system are referred to 
in terms of some new set of coordinate axes. A transformation may be useful in the 
analysis of multiple experiments conducted in different orientations. 

Before any such man ipulation can be performed on data, the data must be collected. 
Various details associated with the sampling interval used in the acquisition of data and the 
window lengths used in the analysis of data are discussed in the next section. The 
adaptation of Fourier theory to the analysis of discrete, aperiodic time series is discussed in 
the Spectral Analysis section. Also discussed in that section are the power spectral density 
function and the application of spectral analysis to residual acceleration data. The 
transformation of data from one set of coordinate axes to another by means of a 
transformation matrix and how this can be useful in the analysis of low-gravity data are 
discussed in the Transformation of Coordinate Axes section. A symbol nomenclature is 
provided in Appendix A, and detailed derivations are included in Appendix B. 


H. SAMPLING CONSIDERATIONS AND FREQUENCY LIMITS 


In the analysis of observational data, certain restrictions are imposed by the length 
of the data window being analyzed and by the sampling rate used when digitizing 
continuous data, or when collecting discrete data. For a segment of a time series f^, of 

length T (N total points), the fundamental period of the segment is assumed to be T, even 
though the series is not necessarily periodic, and the lowest frequency represented in 
spectral analysis of the segment is 1/T = Vi. 1/T also represents the highest resolution 
obtainable in spectral analysis of the time series segment, l/T=Vj=Av. 

The sampling interval At used in the acquisition of data must be appropriate for the 
data-two time intervals are needed to define one period. The sampling interval determines 
the highest obtainable frequency in spectral calculations: v N = l/(2At) = N/(2T), the 
Nyquist frequency. Therefore, for a window of a time series f^ as described above, the 

frequency limits 1/T <,v <, N/(2T) exist for spectral analysis of the window. 

To avoid aliasing, the contamination of computed spectra by frequencies higher 
than v N (see Bendat and Piersol, 1966, pp. 278-280; Waters, 1981, pp. 121-123), 
sampling must be frequent enough (i.e.. At small enough) to have at least two sampling 
intervals (three sample points) per cycle for the highest frequencies present in a series, not 
just the highest frequencies of interest. An investigator must therefore, have some idea of 
the frequency components that will be present in a series before data are collected. Aliasing 
problems may also be avoided by appropriate low-pass filtering as part of the data 
collection process to remove components of frequency higher than the v N dictated by the 
chosen sampling interval. 
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m. SPECTRAL ANALYSIS OF DISCRETE TIME SERIES 
A. Fourier Series and Spectral Analysis - Theory 

Spectral analysis consists of the description of a given time function in terms of 
sinusoidal components present in the function. Such analysis is performed using theory 
developed in part by Joseph Fourier and is referred to as Fourier analysis. According to 
Fourier theory, a function f(t), having a fundamental period of 2 tc and satisfying Dirichlet 
conditions (see Bath, 1974, p. 26), can be represented by an infinite series (a Fourier 
series): 


f(t) = ao + X ( a n cos nt + bn sin nt), 


( 1 ) 


n=l 


where ao, a n , and b n are constants (the Fourier coefficients) which can be represented by: 


Bo = 2k\ *»* 

J-n 

a n = ^-| f(t) cos nt dt, 

n I 

J-n 


(2a) 


(2b) 


and 


-if 


f(t) sin nt dt. 


(2c) 


In general, for a function with fundamental period T, the Fourier series and Fourier 
coefficients are: 


f(t) = ao + X ( a n cos + b n sin 2sni), 


(3a) 


n=l 


ao 


= ^| f(0 dt, 

J -T/2 


(3b) 
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(3c) 


r 

a n = f(t)cos2ffintdt, 

J-T/2 

and 

r 

b n = ^ I f(t) sin dt. (3d) 

J -T/2 

It is the calculation of the Fourier coefficients a n and bn which is referred to as Fourier 
analysis. The value ao is often referred to as the d.c. (direct current) or steady state 
component of f(t), while a n and b n are referred to as the a.c. (alternating current) 
components. 

The Fourier transform (a transformation between the time and frequency domains) 
can be used to calculate the Fourier coefficients of a given time function. The direct Fourier 
transform is derived in the following (derivation in part from B&th, 1974, pp. 35-37, and 
Huang, 1966). The Fourier integral is defined as 


f(t) = 



f(X) cos (co(t - X)) dX, 


(4) 


where X is a dummy integration variable. The cosine and sine transforms of a function f(X) 
arc: 


- 


a(co) = I f(X) cos o )X dX 


(5a) 


and 


■/. 


b(co) = I f(X) sin coX dX. 


(5b) 


Define a function <X>(co) where 


sin <J>(co) = - 


b(co) 


V a 2 (o))+ b 2 (©) ’ 


(6a) 


from which it follows that 


4 



, . . a(co) 

cos O(co) = = == 

V a 2 (co)+ b 2 (co) 


(6b) 


and 


tan O(co) = - 


b(co) 

a(co)' 


(6c) 


Substituting (5a) and (5b) into the Fourier integral (4), with the cosine term 
expanded, results in 
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f(t) = 1 1 [a(co) cos cot + b(co) sin cot] dco; 


(7) 


see Appendix B for derivation. Introducing equations (6) involving O(co) into equation (7) 
yields: 


and also 


f(t) = 1 ( [a 2 (co) + b 2 (co)] 1/2 cos [d>(co) + cot] dco 

n Jo 


f(t) = _L I [ a 2 (co) + b 2 (co)] 1/2 e i<& ( ( °) e iM dco 
2k I 


( 8 ) 


(9) 


= 2nj F(co) eit0t dtD - 


From equations (5), (6), (9), and (10), we can say that 

F(co) = [F(co)| e i<I>(t0) = [a 2 (co) + b 2 (co)] 1/2 e^) 

= [a 2 (co) + b 2 (co)] 1/2 [cos d>(co) + i sin O(co)] 


= [a 2 (co) + b 2 (co)] 1/2 


= a(co) - i b(co) 


a(co) 


b(co) 


[V a 2 (co) + b 2 (co) V a 2 (co) + b 2 (co) 


( 10 ) 


(ID 

( 12 ) 

(13) 

(14) 
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(15) 



f(t) cos tot dt - 1 



f(t) sin cot dt 


f(t) e i<nt dt 


(16) 


From the above, the Fourier transform of a time function (direct Fourier transform, 
Fourier spectrum, complex Fourier transform) is: 


F«o) = 



f(t) e 1(0t dt 


(17) 


where 

F(co) = a(co) - i b(to) = |F(co)| e f *<®>. 

The amplitude spectrum is 

|F(co)| = [a 2 (co) + b 2 (to)] 1/2 


and the phase spectrum is 

O(to) = tan 


-if b(co) 

" a(co) ' 


(18) 


(19) 


( 20 ) 


The Fourier spectrum F(co) of a time function is basically an average of f(t)e ,<Dt over the 
length of the function, i.e., an average of the components of f(t) of frequency co. F(to) and 
f(t) have the same dimensions. 


The real and imaginary parts of F(to), a(co) and b(co), respectively, as shown in (14) 
and (15), are the cosine and sine transforms of f(t) as in equations (5). Comparison of 
these equations to those for the Fourier coefficients, a n and b n (2), shows that the 
following relations exist: 

a(co) an/2 and b(o>) bn/2. (21) 

This can be shown as follows: 
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a(co) = 



f(X) cos coXdX= f f(t) cos 2ma<it 


lim 

T— > oo 



f(t) cos dt = 


The inverse Fourier transform is 


( 22 ) 


f(t) = -L I F(co) e ,(0t dco. 

2k l 


(23) 


This expression can be used to create a time function f(t) from known spectral components 
F(to); therefore, (23) is referred to as the Fourier synthesis of f(t). 

B. Fourier Series and Spectral Analysis - Discretization 

The definitions and derivations given in the previous section are appropriate for 
continuous functions. The calculation of the Fourier transform F(co) of a given function 
f(t) must be approached differently when dealing with a discrete series of observations and 
when processing data on digital computers. The formulas used must be transformed into 
discrete form as must data which are recorded in analog form. For data recorded digitally, 
decisions must be made prior to data acquisition concerning the sampling rate necessary to 
avoid aliasing problems, as discussed in the previous section. See B&th (1974, Chap. 4) 
for more details about the equations presented in this section. 

In general, a continuous integral may be expressed in discrete form as shown here: 


J |*NAx 
0 


N-l N-l 

y(x) dx y(nAx) Ax = ^ y„ Ax, 

n=0 n=0 


(24) 


where y(x) is a given curve, Ax is the digitizing interval, and N is the total number of 
samples. 

An integral with infinite integration limits can be approximated using a form of the 
trapezoidal rule as follows: 
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(25) 


£ 


. X^~X 


J-xn 


N 


y(x) dx -> lim I y(x) dx -* ^ y(nAx) Ax 


NAx n= o 


Assuming y(0) = y(N), the sum of (25) can be represented as: 
iX y( nAx > or y(nAx). 


N 


n=l 


(26) 


n=0 


Using the general forms stated above, the Fourier coefficients and Fourier series 
can be expressed in discrete form as: 


N 

a o = iJ-Z f X> 


(27a) 


a n = ^-S f X cos 
1 


(27b) 


b n = ^-Z4sin2a^, 

Jt=i 


(27c) 


and 


N£ W2 

fx = ao + I a n cos2zffik+ £ bnSin^anA, 

n=l n=l A 

Applying equation (25) to the direct Fourier transform formula yields: 


■L 


N 


F(co) = I f(t) e i0 * dt -> f x c-**. 


(27d) 


(28) 


Equation (28) can be rewritten as follows for ease in computation: 

N 


F n = iS f X e-‘'(2nnX/N). 


X=1 


(29) 
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Comparison of (29) to the continuous Fourier transform formulas (14-16) and the Fourier 
coefficient equations stated above (27) results in: 

Fn=^(a n -ib n ), ( 3 °) 

i.e., the real and imaginary parts of a Fourier transform are one-half the values of the 
Fourier coefficients, an and bn, respectively. The inverse Fourier transform is also easily 
represented in discrete form: 

N 

f X = X F n ei(2nnVN) ' (31) 

n=l 


C. Power Spectral Density Function - Definition 

The power spectral density is a frequency domain function which is often used to 
indicate the dominant frequency components present in data (see Bendat and Piersol, 1966, 
Chap. 3). Before we can define the power spectral density function, we must first 
introduce the autocorrelation function (Robinson and Treitel, 1980, Chapters 3 and 6). The 
autocorrelation function for random data gives some indication of the relative amount of 
power at different frequencies in the data. The autocorrelation function is the expected 
value of the product of a time series and a time shifted version of that series 

R(t) = E(VJ, (32) 

which can also be written as 


R(t) = lim ^ f f(t) f(t+x) dt. (33) 

T->oo 1 ) 0 

The autocorrelation function is a real and positive valued even function with the maximum 
value occurring at T=0. This maximum value corresponds to the mean square value of f(t), 
and represents the energy or power of the time function. 

The power spectral density S(v), associated with the time function f(t), with 
frequencies defined over the interval (- 00 , 00 ), is related to the autocorrelation function: 

R(0) = E(f 2 (t)) = | S(v) dv. (34) 

The power spectral density function represents how the mean square value of a time 
function is distributed over the infinite frequency range. S(v) is called the two-sided power 
spectral density function of f(t) because of the frequency interval (- 00 , 00 ) over which it is 
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defined. It can be shown that the two functions in equation (34) form a Fourier transform 
pair: 


i 


R(x) = S(v) t anvx dv S(v) = R(x) e- i2rcvx dx. 


■I 


(35) 


Because of the nature of the autocorrelation function and properties of the Fourier 
transform, the power spectral density function is a real and positive valued even function of 
frequency. The above transforms may therefore be simplified to: 


■I. 


S(v) = I R(x) cos 2tuvt dx 


-r 

Jo 


R(x) cos 2jtvx dx 


(36) 


and 

r 

R(x) = 2 I S(v) cos 2nvx dv. ( 37 ) 

Jo 

Assuming that R(x) exists, and that it has the Fourier transform S(v), the power 
spectral density function can be defined as the direct Fourier transform of the 
autocorrelation function: 


S(v) = 



R(x) e-' 2KVX dx. 


(38) 


The physically realizable one-sided power spectral density function G(v) for the frequency 
interval (0,°°) is defined as: 

G(v) = 2S(v), 0 < v < oo. (39) 

The correspondence between the one-sided power spectral density function and the 
autocorrelation function is: 


G(v) = 4 



R(x) cos 27tvx dx 


(40) 


and 
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r 

R(x) = I G(v) cos 27CVT dv. 

Jo 


(41) 


D. Power Spectral Density Function - Estimation 

To get the power spectral density function into a form that can be computed directly 
from observational data, without calculation of the autocorrelation function, we must 
consider the estimation of sample parameters (Bendat and Piersol, 1966, Chap. 5). The 
estimation of a previously defined parameter will be indicated by a hat ( A ) symbol over the 
parameter symbol. For example, e is an estimate of e, where e is any parameter, such as 
the mean value or power spectral density of f^, a sample time series existing for the finite 

time interval T. Estimators are often defined arbitrarily and may not clearly give a correct 
estimation of a given parameter. 

Consider the time function fit); the autocorrelation function R(x) can be estimated 
by the sample autocorrelation function R(x), for fit) existing over the time interval T by 

/■ TT 

R(x) = ^J fit) fit+x) dt, OiSxcT. (42) 

Assuming that the data exist for time T+x, the sample autocorrelation function can be 
estimated by 



fit) fit+x) dt. 


0 £ X < T. (43) 


The formula for the estimation of the power spectral density function is not defined 
directly in terms of the autocorrelation function, but takes into account the fact that the 
autocorrelation function is related to the mean square value of a time function. For a time 
function fit) with zero mean, existing over a time interval T, the estimate of the one-sided 
power spectral density function G(v) describes the time average of f 2 (t) in terms of its 
frequency components, in the frequency interval (v-(Be/2), v+(Be/2)), where B e is the 
frequency bandwidth. 

The mean square value of fit) within the bandwidth B e centered at v is estimated by 


? 2 (v3e) = 



f ^t.v.Be) dt. 


(44) 
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The power spectral density is defined as 


G(v) = lim 

B,— »0 


y 2 (V,B e ) 

B c 


lim 

T— » oo 
B e — >0 



f 2 (t,V,Be) dt 


= lim G(v). 

T— » «• 

B e -»0 


( 45 ) 


From equations (44) and (45), the sample estimate for the one-sided power spectral density 
function G(v) is: 


d(V >-' L W A ^[ ( 46 ) 

The limits in equation (45) must be evaluated to obtain the true function G(v); G(v) is 

/s 

generally a biased estimate of G(v). The power spectral density function estimate G(v) is 
one-sided and is related to the mathematical two-sided power spectral density function 
estimate defined for positive and negative frequencies as shown here: 

S(v) = S(-v) = (47) 

The above formula (46) for the estimate of the power spectral density function can 
be discretized as follows, taking into account that the bandwidth B e is equivalent to the 
Nyquist frequency v^: 



_ 1 
VnN 



x=o 


1 


N-l 


(N/2T)N 


I* 

x=o 


(48a) 


(48b) 


(48c) 
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(49) 


N-l 


= 2T V 

N 2 ^ 
^ X=0 


f 2 

Iv • 


The estimate of the one-sided power spectral density function as defined above is 
related to the Fourier spectrum of a time series. Working from the formula for the discrete 
Fourier transform of f^, F n =(a n -/bn)/2, where a n and bn are the Fourier coefficients of f^, 

the relation can be shown as follows: 

F n = 2 ( a n - * bn) ( 5 °) 



N 2 


N N 

SI 

U=l X=l 


fx 2 COS 2 



sin z 




X X ^ 2 ( c °s 2 ( 2 ^ ) + Sin 2 ( 27 J^) ) 
Jl=i X=i 


(51a) 

(51b) 

(51c) 

(5 Id) 
(51e) 


N 




(52) 


Comparing equation (52) with the discrete equation for the estimate of the power spectral 
density (49) yields the following relation: 

G n m 2TF 2 . ( 53 ) 

This can be rewritten, to ease comparison to the Fourier transform of a time series of 
interest: 


13 


2 


PSD n = G n = 2E- 
N 2 


N 


y' e -i2jmVN 
Jl=l 


( 54 ) 


The units of the power spectral density function are (units of Fourier transform) 2 /Hz. 

E. Application to Residual Acceleration Data 

Residual acceleration data have been collected in orbiting space laboratories with a 
variety of instruments. The data are typically recorded and distributed to investigators in 
discrete form. Sampling rates vary considerably from one experiment to the next, 
depending on the specific goal of the experiment. Typical sampling rates range from 12.5 
to 500 samples per second. Sampling is often done at higher rates than required to obtain a 
particular maximum frequency. For example, if the highest frequency of interest to an 
investigator is 100 Hz, data may be collected at 500 samples per second and then lowpass 
filtered down to 100 Hz. Such a procedure is useful when investigators are not sure of the 
maximum frequencies that compose the process being sampled. This also allows for 
higher resolution of frequency domain data than would be available if a sampling rate of 
only 200 Hz was used. 

Before any specific Fourier analysis can be done on accelerometer data, some pre- 
processing is usually required. Pre-processing is performed so that the output represents 
as closely as possible the actual residual accelerations experienced at the recording site. 
Artificial signals, or noise, can be introduced from a multitude of sources. Manufacturers 
of recording devices typically furnish users with specific corrections for deviations from 
pure signal caused by temperature variations, instrumental bias, and other factors. Filtering 
can also be applied to data to minimize the effects of aliasing, instrumental noise, and other 
known noise sources. Specific filtering techniques wifi not be discussed here, but 
discussions of filtering can be found in most signal processing texts, see B&th (1974), 
Bendat and Piersol (1980), and Cadzow (1987). 

Once such pre-processing has been applied to residual acceleration data, Fourier 
analysis can be applied in order to determine what frequency components are significantly 
present. Windows must be chosen of a length appropriate to what minimum frequencies 
one wants to see, as discussed in the Sampling Considerations section. Many Fourier 
transform algorithms are most efficient when the number of samples is highly composite, 
e.g., a power of two, which places further restrictions on the window length. After pre- 
processing, the data are assumed free of artificial signals. The results of inadequate 
correction for instrumental bias can often be seen in Fourier transformed data. An artificial 
bias from the zero point will cause a significant steady state component. Other incorrect 
processing, however, is not as easy to identify, but may be manifested as a simple bias, or 
as long-term variations or trends. 

Most standard Fourier transform computer programs use algorithms based on the 
discrete Fourier transform formula 
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( 29 ) 



and may give as output the cosine and sine transforms, a(co) and b(co), respectively, or the 
amplitude and phase spectra: 


F n = V a 2 (co) + b 2 (co) 


— tan’l 


The Fourier coefficients a n and bn which differ from the cosine and sine transforms by a 
factor of 2 may also be output, so it is important to completely understand the transform 
routine being used. See, for example, Bloomfield (1976, Chap. 4), Elliott and Rao 
(1982), and Cadzow (1987, Chap. 4) for information on Fourier transform algorithms. 

Given the Fourier transform F n , the original series may be recreated using the 
discrete inverse Fourier transform (31) or it may be synthesized using the relation between 
F„ and the Fourier coefficients (30) and the discrete Fourier series equation (27d). Using 
equation (54), output from Fourier transform programs can be used to form the power 
spectral density function of the time series considered which can be investigated in addition 
to, or as an alternative to, the Fourier spectrum. 
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IV. TRANSFORMATION OF COORDINATE AXES 


A. Transformation of Coordinate Axes - Theory 

Observational time series data are usually recorded in a stationary set of coordinate 
axes, which may be located and oriented arbitrarily. Such data can be analyzed in terms of 
alternative axes by a relatively simple transformation operation, as long as the 
transformation parameters between the axes of interest and the recording axes are known. 

One common form of coordinate axes transformation is that which transforms axes 
that share a common origin and differ by an angle of rotation 0 about a common axis; see 
Figure 1. In the transformation of axes, a known vector g, such as a residual acceleration 
measurement, is written in terms of some (primed) coordinate axes which are related to the 
original (unprimed) coordinate axes by a known rotation angle. That is, g has some 
coordinates in the unprimed axes (gi,g2,g3). and some other coordinates in the new, 

primed axes (gi'^'.gaO- 



Figure 1. Axes a .2 and are rotated by an angle of 0 about the ai axis, g is some vector 
with coordinates (gi, g 2 , g3) in the unprimed axes and (gi', g 2 ', g 3 7 ) in the primed axes. 
Unit vectors e and e' are marked by short vectors along the coordinate axes. 
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Referring to the vector g in terms of some new coordinate axes is equivalent to a 
change in basis of a vector space from an original basis to a new basis (Anton, 1981, 
Chap. 4). The vector g is rewritten as g' in terms of the new coordinate axes, g' is 
obtained by the application of a transition matrix to g: 


g' = Rg 


(56) 


R is called the transformation matrix (or rotation matrix). The columns of R are the 
components of the old basis vectors relative to the new basis. Consider, for example, the 
unit vectors ei, *i, ej\ and of two coordinate axes rotated with respect to one another 
by an angle 0. The transition matrix can be constructed by writing ej and e 2 in terms of 
the primed axes: 


ei = ei' cos 0 - e 2 ' sin 0 
e 2 = ei' sin 0 + cos 0 


which yields the transition matrix: 


R = 


cos 0 
-sin 0 


sin 0 

cos 0 . 


This is the transition matrix from unprimed to primed coordinate axes: 


g' = Rg. 


(57) 


(58) 


(59) 


This can be checked by writing out the elements of g' : 


gi' = gi cos 0 + g 2 sin 0 
Z2 = -gi sin 0 + g 2 cos 0 


(60) 


and comparing them to the equations for ex' and in terms of the unprimed axes: 


ex' = ex cos 0 + e 2 sin 0 
- -ei sin 0 + e 2 cos 0. 


(61) 


Equations (60) and (61) are equivalent. Equations (61) could be used to create the 
transition matrix for primed to unprimed coordinate axes, which would be equal to the 
inverse (R- 1= R0 of the transition matrix from unprimed to primed coordinate axes. The 
process is easily extended to three dimensions (Anton, 1981, Chap. 4). 
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In general, for a coordinate system defined by three mutually perpendicular axes, 
the transition matrix can be constructed from the direction cosines of the primed axes with 
respect to the unprimed axes: 



Rn 

Rj 2 

Rl 3 

R = 

R21 

R22 

R23 


-R31 

R32 

R33- 



(62) 


where the first subscript denotes a primed axis and the second an unprimed axis (Frederick 
and Chang, 1972, Chap. 1). For example, R 12 is the cosine of the angle between the ai'- 
axis and the a 2 -axis. Hence, the transition of a vector from an unprimed system to a 
primed system can be obtained by 

gi' = Rij gj- ( 63 ) 

The use of equation (62) is preferred over that of equation (58) because (62) is more 
general and can be used for cases where the rotation is around any line through the origin, 
not only around a shared axis. 

B. Application to Residual Acceleration Data 

The transformation of coordinate axes can be a very useful operation in the analysis 
of data. In the case of the analysis of residual acceleration data in conjunction with 
experimental data, the ability to consider data in terms of different coordinate axes is a 
necessity. Experiments run in orbiting space laboratories, including accelerometer 
systems, are often oriented in a manner convenient to space restrictions. This results in a 
plethora of experimental coordinate systems for one mission. 

While the total acceleration vector has the same magnitude in any coordinate 
system, the disturbance level in some particular direction may be of interest to 
investigators. Some experiments may be more sensitive to disturbances in one direction 
than in another direction. In the analysis of the results of such experiments, it is beneficial 
to transform the accelerometer data into a set of axes coincident with the direction of 
interest. The magnitude of disturbances in the direction of interest can then be obtained. If 
the sources of disruptive disturbances can be identified, and if these disturbances tend to be 
uni-directional, such knowledge could be used in the future when positioning experiments 
in orbiters. 
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V. CONCLUSIONS 


Various aspects of residual acceleration data are of interest to investigators. The 
mean, variance, and minimum and maximum values can be obtained from the data as 
collected in the time domain. Additional information can be obtained by looking at data 
from a different perspective. Dominant frequency components can be identified following 
a transformation of data into the frequency domain (amplitude spectrum or power spectral 
density) using Fourier transform methods. Information obtained in both the time and 
frequency domains can be used by investigators to determine what magnitude disturbances 
and what frequency modes are disruptive to their experiments. 

The orientation of accelerometer recording axes and of recorded accelerations is 
important in the analysis of low-gravity experiments. The ability to refer to acceleration 
data in terms of different coordinate axes (such as those of a separate experiment) is useful 
in the post-flight analysis of experiments. 

Methods such as those discussed here can be implemented in the analysis of 
residual acceleration data collected in orbiting space laboratories and used to support 
analysis of experiments run under low-gravity conditions. 
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APPENDIX A - NOMENCLATURE 


a„, bn Fourier coefficients 

a(o)) cosine transform of time function 


b(co) sine transform of time function 

ei, e2 unit vectors 

f(t) function of time, t 

f^ time series with increment K, digitized version of f(t) 

F(©) Fourier transform of f(t) 

F n Fourier transform of f^ 

g vector with coordinates (gi, g2, g3) 

G(v) one-sided power spectral density function 

i VT 

N length of time series in number of samples 

Rij transformation matrix = R 

R(x) autocorrelation function of a time series 
S(v) two-sided power spectral density function 
T length of time series or function in seconds, 

also fundamental period of time series or function 
At sampling interval 

$(©) phase spectrum of f(t) 

0 angle of rotation between coordinate axes 

x time shift 

v cyclic frequency 

v N Nyquist frequency 

© angular frequency, ©=2jcv 

A over a symbol in the text represents an estimate of the parameter. 
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APPENDIX B - DET AILED DERIVATIONS 


The derivation of equations (9) and (10) in the main text from the Fourier integral 
(4) can be done as follows. Given the Fourier integral 


f(t) 


= *I **! 


f(X) cos (to(t - X)) dX, 


(Bl) 


the cosine term can be expanded, yielding: 


f(t) 


■HH 


f(X) [cos ©t cos coX + sin tot sin coX] dX 


(B2) 




[f(X) cos cot cos toX + f(X) sin tot sin toX] dX 


(B3) 


i: 


f(X) cos coX cos tot dX + I f(X) sin toX sin tot dX 


(B4) 


-if 

’'Jo 


[a(to) cos tot + b(to) sin tot] dto. 


(B5) 


Equation (B5) is the same as equation (7) in the main text. Introducing equations (6) from 
the main text for d>(to), a(to), and b(to) yields: 

f(t) = -X I [a 2 (to) + b 2 (to)] 1/2 cos O(to) cos tot 

Mo 

- [a 2 (co) + b 2 (to)] I/z sin O(to) sintot dto (B6) 


-if 

"Jo 


[a 2 (co) + b 2 (to)] 1/2 [cos O(to) cos tot - sin O(to) sin tot] dto (B7) 
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-if 

”1 


[a 2 (co) + b 2 (co)] w [ COS (*(G» + cot)] dco. 


(B8) 


Equation (B8), which is equation (8) in the main text, can he rewritten as follows: 

f(t) = -1- I [ a 2(co) + b 2 (co)] 1/2 [eW ffl > +®«> + e-»W®> + ® t )] dco 

27t io 


(B9) 


■if 


[a 2 (co) + b 2 (co)] 1/2 [eW®> +<Dt >] dco 


27C 


/•o 

I [a 2 (co) + b 2 (co)] 1/2 [e‘‘W®> + «>] dco. 


(BIO) 


by reversal of the limits of integration. Manipulation of the exponential terms yields: 


f(t) = -L I [ a 2 (co) + b 2 (co)] 1/2 [e*® 1 e^H dco 

2n Jo 


f° 

-L [a 2 (co) + b 2 (co)] 1/2 [e‘'«* e*'^®)] dco. 

2n L 


(fill) 


Addition of the two integrals results in equations (9) and (10) of the main text, the inverse 
Fourier transform: 


/ °° 


f(t) = ^j [a 2 (co) + b 2 (c o)] 1/2 [e 1 ® 1 eW«>] dco 


(B12) 


’^i 


F(co) e*® 1 dco. 


(B13) 
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